QND measurement of a superconducting qubit in the weakly projective regime 
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Quantum state detectors based on switching of hysteretic Josephson junctions biased close to 
their critical current are simple to use but have strong back-action. We show that the back-action 
of a DC-switching detector can be considerably reduced by limiting the switching voltage and using 
a fast cryogenic amplifier, such that a single readout can be completed within 25 ns at a repetition 
rate of 1 MHz without loss of contrast. Based on a sequence of two successive readouts we show 
that the measurement has a clear quantum non-demolition character, with a QND fidelity of 75 %. 
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According to the Heisenberg uncertainty principle, the 
measurement of a physical variable necessarily perturbs 
its quantum conjugate variable, imposing a fundamental 
limit on the precision of the measurement. However, it 
is in principle possible to decouple the dynamics of the 
variable which is measured from the perturbed conju- 
gate variable. This is the central idea in the concept of 
Quantum Non-Demolition (QND) measurement, which 
was initially developed to perform high precision mea- 
surement beyond the standard quantum limit [1 . QND 
measurement became an important paradigm of quantum 
mechanics of particular interest for fundamental studies 
of quantum systems. A simple criterion ensuring that 
a measurement is QND is that the Hamiltonians of the 
bare quantum system and of the interaction between the 
quantum system and the detector do commute. In prac- 
tice, its realization is not straightforward as it requires a 
very good control of the back-action of the classical de- 
tector on the state of the quantum system. So far, QND 
measurements have been demonstrated experimentally in 
quantum optics [2], in cavity-QED with Rydberg atoms 
[3], with trapped single electrons [4^, and more recently 
with superconducting circuits [5HZ] based on dispersive 
coupling of a quantum 2-level system (or qubit) with a 
resonator- type detector. 

In this letter we study the quantum measurement of 
a superconducting flux qubit, employing a hysteretic dc 
SQUID, which is either in a superconducting state or in 
a dissipative state. The back-action of this detector, in 
particular if in the dissipative state, is commonly much 
higher than the minimum imposed by quantum mechan- 
ics, thus seriously limiting its efficiency. In our experi- 
ment, we show that a single readout can be completed in 
a time as short as 25 ns, and that the back-action can be 
controllably reduced to such a level that for the first time 
with this type of detector, a QND measurement is pos- 
sible with an averaged QND fidelity of 75 %. It should 
be stressed that although this type of fast DC-switching 
detector can generally provide high readout contrast [8|, 
in this experiment we chose a small qubit-detector inter- 
action resulting in a low readout contrast and therefore 
a weakly projective measurement. This simple and fast 
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FIG. 1: (color on-line) (a) Scanning electron microscopy pic- 
ture of the flux qubit and the dc SQUID, (b) Optical picture 
of the readout circuit, (c) Schematic of the readout circuit 
including a dc SQUID shunted with a resistor R — 125 Q 
and a capacitor C — 1 pF, and a HEMT cryogenic amplifier 
allowing a very fast detection of the qubit state. 



QND detector forms a highly attractive alternative for 
dispersive QND detectors. 

The measurement can be simply illustrated in the 
Bloch sphere representation of the qubit state. Denot- 
ing as 1^ > and |e > the ground state and excited state 
respectively, the general qubit state can be written as 
\ip >= cos (f)l^ > +6*"^ sin (f)|e >. The variable which 
is measured is the qubit excited state occupation prob- 
ability p{e) = sin^ |, while ^ is the perturbed conjugate 
variable. As the free evolution of 6 is decoupled from 0, 
the main requirement for a QND measurement is that no 
transitions between \g > and |e > occur during the mea- 
surement. Due to the small discrimination efficiency be- 
tween |e > and |^ > in a single readout, p{e) is obtained 
from the averaged outcome of a large ensemble of iden- 
tical qubit state preparation and readout sequences. In 
the case of an ideal QND detector, the averaged outcome 
of two successive measurements should be equal. Ad- 



ditionally, the correlation of the measurement outcomes 
reveals the degree of projectiveness of the measurement 
which in our experiment is very small. 

The flux qubit, shown in Fig. [l^, is made of a super- 
conducting loop interrupted by four Josephson junctions 
(the largest junction being auxiliary), biased with an ex- 
ternal flux ^ close to half a flux quantum ^*^o = h/Ae [9]. 
In general, the two energy eigenstates \g > and |e > are 
quantum superpositions of two states corresponding to 
oppositely circulating persistent currents of magnitude Ip 
and denoted as |+ > and |— >. In the basis of these cur- 
rent states the Hamiltonian of the flux qubit can be writ- 
ten as: H = —^{edz + A(Xe), where <Xe, dz are Pauli spin 
matrices, e = 2/p($ — <l>o/2) and A is the tunnel coupling 
between the two current states. The difference of energy 
between \g > and \e > is E = Ee - Eg = -J e^ + A^. All 
the measurements are performed at a bias |e| >> A, 
where the energy eigenstates \g > and |e > approxi- 
mately coincide with the current states |+ > and |— >. 

The quantum state of the qubit is read out by measur- 
ing the flux resulting from the expectation value of the 
circulating current, using a hysteretic dc SQUID compris- 
ing two Josephson junctions in a superconducting loop. 
The interaction between the qubit and the detector is 
given by Hj = MIp&zJ^ where M is the qubit-dc SQUID 
mutual inductance and J is the circulating current in the 
SQUID loop. We should stress that the mentioned simi- 
larity between the eigenstates of the current and the en- 
ergy for \€\ » A closely approaches the sufficient condi- 
tion for a QND measurement [H, Hj] = 0. Furthermore, 
the large difference between the internal frequency of the 
detector and the qubit frequency suppresses energy ex- 
change between qubit and detector. 

The dc SQUID is operated as a DC-switching detec- 
tor which is either in a purely superconducting state S 
with zero voltage, or in a dissipative state D with a finite 
voltage. The dc SQUID, initially in the superconducting 
state S is biased with a current I^ close to its critical cur- 
rent Ic with a short current pulse where it might switch 
to the dissipative state D with a probability P. The state 
of the dc SQUID is maintained using a lower bias current 
h ^ ^c/2, for a time allowing to discriminate efficiently 
the two dc SQUID states. Figures |2^ and ^p show the 
readout current pulse sent to the dc SQUID and the re- 
sulting output voltage in the case where the dc SQUID 
is in the state S (cyan) or switches to the state D (ma- 
genta). P depends on the flux through the dc SQUID 
loop and therefore on the state of the qubit. The de- 
tector switching probabilities in the two cases when the 
qubit is in the ground state or in the excited state are 
denoted as P{g) and P{e). For a general qubit state with 
an excited state occupation probability p{e) and a read- 
out contrast Cm = P{^) — P{9) the switching probability 
is 
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FIG. 2: (color on-line) (a) Current pulse for the qubit read- 
out, (b) Output voltage. Dots correspond to real time signal, 
while the solid lines are average of 16 traces corresponding to 
the state S or D. (c) Histogram of the output voltage, (d) 
Energy spectrum of the flux qubit. The color scale indicates 
the switching probability of the detector, (e) Rabi oscillation 
for a measurement repetition rate of 10 kHz or 1 MHz with 
an identical manipulation and readout sequence at a qubit 
frequency of E/h=13 GHz. The three dashed lines represent 
P(g), P(e) and P^. 



Switching of the dc SQUID to its dissipative state has 
several severe consequences. The first trivial but impor- 
tant effect is that Joule heating in the shunt resistor 
increases the temperature of that resistor, which leads 
to additional noise. The dissipated energy also leads to 
an increased temperature of substrate and qubit. When 
the dc SQUID is at a nonzero voltage, the AC Joseph- 
son effect gives rise to currents at frequencies of tens to 
hundreds of GHz, which couple into the qubit and may 
cause non-adiabatic transitions and quantum leakage to 
higher levels. If the voltage is above the superconducting 
gap, generation of quasiparticles occurs in the dc SQUID 
junctions, deteriorating the superconducting properties. 
Quasiparticle generation in the qubit by transferred RF 
fields, and by direct diffusion when the qubit is gal- 
vanically attached, is detrimental to qubit performance. 
Quasiparticle recombination takes a long time and very 
low numbers of quasiparticles in a qubit can be signifi- 
cant. All these effects lead to long recuperation times for 
dc SQUID and qubit between successive readouts, com- 
monly limiting the measurement repetition rate to a few 
tens of kHz. 

To reduce the back-action of the detector, we use a 
very simple approach [lOl |TT] which has never been in- 
vestigated so far for non-destructive readout of supercon- 
ducting qubits. By shunting the dc SQUID with a low 
external resistor R^ the voltage in the dissipative state 
can be reduced to a value of V ~ RI^ much below the 



superconducting gap. Therefore the generation of heat 
and quasiparticles is greatly suppressed. Heating is fur- 
ther reduced by using a fast cryogenic amplifier allowing 
to discriminate the two dc SQUID states efficiently with 
a very short readout pulse of 25 ns duration. To avoid 
non-adiabatic transitions of the qubit caused by high fre- 
quency oscillation of the current in the dc SQUID due 
to the AC-Josephson effect, the critical current of the dc 
SQUID is chosen such that the qubit-detector interaction 
is smaller than 1 GHz. 

The sample is fabricated using standard e-beam lithog- 
raphy and shadow evaporation (Fig. Ilk andfTb). On- 
chip film resistors are made of AuPd while all supercon- 
ducting structures are fabricated of aluminium, with a 
thin aluminium oxide layer to form the Josephson junc- 
tions. The dc SQUID has two junctions of critical cur- 
rent Ico = 80 nA and is shunted by an interdigitated 
finger capacitor C = 1 pF and a resistor R = 125 Q. 
It is operated at a critical current Ic = 140 nA, lead- 
ing to F ^ RIc = 17 /iV. The qubit is characterised by 
Ip = 400 nA and A = 3.5 GHz. It is directly coupled 
to the dc SQUID via a shared narrow line section with 
a mutual kinetic inductance M = 30 pH. The maximum 
interaction between the qubit and the dc SQUID during 
the measurement is 2MIpJ/h = 0.7 GHz. 

The dc SQUID output voltage is measured by an ultra- 
low noise cryogenic HEMT amplifier (Fig. flfc), charac- 
terised by a high input impedance, 0.3 GHz bandwidth 
and ^ 0.25 uV/a/Hz noise. The amplifier is mounted on 
a Cu strip attached to the refrigerator still at ~ 600 mK. 
Its input is connected to the sample by two 7 cm length 
superconducting wires. Rout = 3 kQ and Cqut — 20 
pF define an overall measurement response time of 60 
ns. Based on the bandwidth, the dc SQUID voltage am- 
plitude and the amplifier noise, the two states of the dc 
SQUID can be discriminated to better than 97 % confi- 
dence with a readout pulse of 25 ns total duration (Fig. 
^). To gain sufficient statistics typically 25000 indi- 
vidual events are taken for each measurement condition. 
The experiments were performed in a cryogen-free dilu- 
tion refrigerator, at a base temperature of 20 mK. 

The resistor R shunting the dc SQUID induces relax- 
ation in the qubit. However, during the coherent ma- 
nipulation of the qubit, this relaxation channel is closed 
by operating the dc SQUID at a bias current such that 
jYb(-^b) = 0, where external ffuctuations of I^ do not cou- 
ple to the dc SQUID circulating current J [TT1[T2]. For a 
symmetric dc SQUID, this condition is realized at /& = 0. 
During the readout, jjiih) 7^ and the qubit relax- 
ation rate when the dc SQUID is in its superconducting 
state S can be estimated using Fermi's golden rule as- 
suming a Johnson-Nyquist spectral noise density in the 
high-frequency limit hu > ksT 
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FIG. 3: (color on-line) (a) Qubit operation and readout se- 
quence including two Rabi oscillations of angles Oa and Ob^ 
and two readout pulses A and B. The qubit operation fre- 
quency is set at E/h=13 GHz. (b) Switching probabilities 
Pa , Pb in readout A and B as a function of the delay be- 
tween the two measurements. P^^^^ is the free qubit relax- 
ation decay measured with readout B without applying the 
readout pulse A. (c) Pa,Pb as a function of the duration of 
the first Rabi pulse Tra- Pb{t^) and Ps(0) are the results 
readout B with or without a tt rotation between the two mea- 
surements. (d,e,f) Switching probabilities in the first readout 
Pa and in second readout in the case that the first readout 
did not Pb{Sa) or did switch Pb{Da) as a function of the 
durations Tra and Trb of both Rabi pulses. 



Based on ([2]), the relaxation rate due to the shunt resis- 
tor during the holding plateau is given by l/F^ c^ 500 ns, 
while during the switching plateau l/F^ ::^ 40 ns. Given 
the total duration of the readout (25 ns) and the width 
of the switching plateau (~ 2 ns), the relaxation dur- 
ing the readout should be small when the detector is in 
state S. When the dc SQUID is in its dissipative state 
D^ qubit relaxation is more complex to estimate as it 
involves strongly non-linear effects. 

We first investigate the back-action of the resistively 
shunted dc SQUID on the coherence of the fiux qubit by 
performing coherent Rabi oscillations at different mea- 
surement repetition rates. As can be seen in Fig. |2^, 
virtually no reduction of coherence is visible up to a mea- 
surement repetition rate of 1 MHz, nearly two orders 
of magnitude higher than previously achieved with un- 
shunted or weakly shunted DC-switching detectors [13]. 
Although the Rabi coherence time of this qubit sample is 
relatively short, the fact that such a high measurement 
repetition rate is possible without any loss of readout 
contrast and qubit coherence shows that the back-action 
is very low. 

The contrast of the readout Cm = P{^) — P{9) is ob- 
tained from the Rabi oscillations in Fig. [2J3. As the qubit 
energy E /h = 13 GHz is much higher than the thermal 



energy kT/h ^ IGHz, the qubit is initially in the ground 
state, corresponding to a detector switching probability 
P{g). We assume that the oscillation of P as a function of 
the duration of the Rabi pulse Tmw is symmetric around 
Pm = {Pg + ^e)/2, such that in the limit of a long and 
incoherent pulse P — Pm- Based on the values of P{g) 
and Pm extracted from Fig. [2J3, we deduce Cm = 12%. 

We now turn to the analysis of the QND character 
of the measurement using a sequence of two successive 
readouts A and B and two Rabi pulses (Fig. Isk). First, 
the switching probabilities Pa and Pb in the two read- 
outs are measured as a function of the delay between the 
two measurements when the qubit is initially prepared in 
the excited state using a Rabi pulse Ba = tt. As shown 
in Fig. [3J3, Pa is constant while Pb decays exponen- 
tially. To quantify the effect of readout A on the qubit, 
we also measure PJ^^^ in the same conditions as previ- 
ously except that the readout pulse A is not applied. The 
exponential decay of Pb^^i which is only due to the free 
relaxation of the qubit Ti=160 ns can be directly com- 
pared with the decay of Pb- As can be seen in Fig. ^p 
the decay time of Pb is equal to the free relaxation time 
(solid lines are exponential decay fitting curves). From 
the ratio Pb/Pb^^ = 0.8 we conclude that the qubit re- 
laxation is increased by an amount of 20 % during the 
measurement due to the operation of the detector. 

In addition to the enhancement of the qubit relaxation 
during the readout, we observe that the result of the first 
readout influences the result of the second readout for 
delay between the measurements shorter than 20 ns. For 
a delay of 10 ns, the switching probability Pb{Da) in 
readout B if the detector was in state D during readout 
A is increased by about 10 % compared to Pb{Sa) if 
the detector was in state S. As this effect is not related 
to the qubit, it can be easily corrected for by measuring 
the difference Pb{Da) — Pb{Sa) with the qubit in one 
specific state, e.g. the ground state. In Fig. [3J3 the open 
dots represent Pb after applying this small correction. 

Secondly, we measure Pa and Pb as a function of the 
duration Tra of the Rabi rotation pulse preceding read- 
out A, with a second Rabi pulse between readout A and 
B, and a fixed delay time of 100 ns between the two mea- 
surements (see Fig. Isk). Pb{t^) and Pb(0) are obtained 
either with or without a tt pulse before readout B. Rabi 
oscillations are visible in Pa, Pb(0) and Pb{^) in Fig. pp. 
The reduced oscillation amplitude of Pb (0) compared to 
Pa is in good agreement with the free qubit relaxation 
during a time of 100 ns and the extra relaxation due to 
readout A mentioned above. The fact that the sign of 
the oscillation of Pb{7t) is reversed compared to Pb(0) is 
a clear proof that the qubit state is partially preserved 
and can be manipulated coherently after a measurement. 
Furthermore, the averaged value of Pb{7t) is higher than 
Pe(0). This is to be expected as the effect of the tt pulse 
is to reverse the qubit state and the qubit relaxation be- 
tween the measurements is significant. 



Figures [3]i,|3^ and[3]F show Pa and the two conditional 
switching probabilities Pb{Sa) and Pb{Da) in readout 
B in the cases that the result of readout A is S or D 
respectively, as a function of Tjia and Tjib- To improve 
the visibility of the oscillations of Pb{Sa) and Pb{Da) 
a constant offset is subtracted for each duration of the 
second Rabi pulse Trb- As the discrimination efficiency 
between \g > and |e > in a single readout is low, the qubit 
state is only very weakly projected after a measurement. 
Given the readout contrast Cm = 12 %, the correlation 
between readout A and B due to the qubit state projec- 
tion is too small to be resolved. The two readouts are 
thus nearly independent. As can be seen in Fig. |3^ and 
[Sf, the qubit state is almost as well preserved if the de- 
tector is in state S or in state D. This result shows that 
the back-action of this type of detector in the dissipative 
state can be reduced to a very low level. 

We define the QND fidelity based on the conditional 
probabilities that the qubit state after the measurement 
is the same as before the measurement, irrespective of 
the measurement outcome, when the qubit is initially in 
1^ > or |e >, denoted as p{g\g) and p{e\e) respectively 
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This definition, which does not involve the readout con- 
trast, is only relevant if the readout contrast is high 
enough compared to the available statistic. Using rela- 
tion (IT]) for the readout B, equation ([3| can be rewritten 
as Fqnd = ^ + ^^^ where APb = PB{e)-PB{g) is the 
difference of switching probability in readout B depend- 
ing on the initial qubit state before readout A. From the 
measurement of Fig. ^p for a delay between the measure- 
ment of 10 ns we find APb > 6 %, without correcting for 
the initial qubit preparation errors but subtracting for 
the systematic extra switching probability due the read- 
out A as explained previously. Assuming the contrast 
Cm = 12 % for readout A as well as for readout B we 
conclude that Fqnd ^ 75 %, with p{g\g) — 100 % and 
p{e\e) > 50 %. 

We should emphasize that this letter presents the re- 
sults of the first attempt to reduce the back-action of the 
DC-switching readout. The readout contrast can still be 
largely improved while maintaining a good QND fidelity. 

In conclusion, we have demonstrated that the total 
readout time and the back-action of the DC-switching 
detector can be considerably reduced, such that a sin- 
gle readout of a superconducting flux qubit can be com- 
pleted within 25 ns with an averaged QND fidelity of 
75 %. Such simple, fast and low back-action detector 
is very attractive to study entanglement in a multiple- 
qubits and multiple-detectors system and opens new per- 
spectives to investigate fundamental aspects of quantum 
measurement. 
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